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Recent phenomenological work has examined two different ways of including charge symmetry 
violation in parton distribution functions. First, a global phenomenological fit to high energy data 
has included charge symmetry breaking terms, leading to limits on the magnitude of parton charge 
symmetry breaking. In a second approach, two groups have included the coupling of partons to 
photons in the QCD evolution equations. One possible experiment that could search for isospin 
violation in parton distributions is a measurement of the asymmetry in W production at a collider. 
In this work we include both of the postulated sources of parton charge symmetry violation. We show 
that, given charge symmetry violation of a magnitude consistent with existing high energy data, the 
expected W production asymmetries would be quite small, generally less than one percent. 



Charge symmetry represents a specific form of isospin 
invariance (a rotation of 180° about the "2" axis in 
isospin space) that is quite well respected at low ener- 
gies pi • Since there is no direct experimental evidence 
of charge symmetry violation (CSV) for PDFs 0,0, it 
was reasonable, at least in the beginning, to assume that 
it held as well for parton distribution functions (PDFs). 
However, we know that small violations of charge sym- 
metry do arise from both the mass differences of light 
current quarks, and from electromagnetic effects. There 
have been some theoretical estimates of charge symmetry 
violation in PDFs, and recently charge symmetry viola- 
tion has been included in phenomenological PDFs. Fur- 
thermore, the estimated size of the CSV is such that it 
can produce important effects in some experiments, for 
example in precise tests of physics beyond the Standard 
Model [10. 

Global fits of PDFs by Martin, Roberts, Stirling and 
Thorne (MRST) included the possibility of charge 
symmetry violating PDFs for valence and sea quarks. 
By construction, the resulting parton distributions will 
agree with the array of experimental data used in global 
fits. The valence quark CSV PDFs were chosen to have 
the specific form 



Su v (x) = v%(x) — d™(x) ; 6d v (x) = d p ,(x) — u"(x) 



(1) 



At both small and large x the valence quark CSV term is 
qualitatively similar to phenomenological valence quark 
distributions , and the first moment of the valence CSV 
distribution is zero, a necessary condition to preserve va- 
lence quark normalization. The single coefficient k was 
varied in the global fit to high energy data. To mini- 
mize the resulting computing time, MRST neglected the 
Q 2 dependence of this CSV effect. The best value they 
obtained was k — —0.2, and the 90% confidence level 
included the range —0.8 < k < +0.65. It is interest- 
ing to note that for the best fit obtained by MRST, the 
valence quark CSV distributions are in very good agree- 
ment both in sign and magnitude with predictions from 
quark model calculations j^, [Hi ~ see a l so the- model in- 
dependent constraint on the second moment obtained in 
Refs. PHG2. 

In a separate global fit to the same data, MRST in- 
cluded the possibility of sea quark CSV effects. The 
MRST functional form chosen for sea quark CSV was 



6u v (x) 



-Sdy(x) = «(1 - xfx-''-" (x - .0909) 



u n {x) = d p (x)[l + S] 
d n (x) = u p {x)[l-S] 



(2) 
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This form was chosen to insure that the total momen- 
tum carried by antiquarks in the neutron and proton was 
approximately equal. Once again, they assumed no Q 2 
dependence for these CSV distributions. The best fit was 
obtained for S = 0.08. 
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An alternative, phenomenological approach to the 
problem of charge symmetry violation associated with 
the electromagnetic interaction has been proposed by 
both MRST UJ and Glueck, Jimcnez-Delgado and Reya 
[l4| . By analogy with the usual QCD evolution involving 
gluon radiation, these authors suggested that one assume 
charge symmetry at some initial low-mass scale, and in- 
clude in the evolution equations the effect of photon ra- 
diation. When one includes QED contributions in this 
way, to lowest order in both the strong coupling a s and 
the EM coupling a, the so-called DGLAP evolution equa- 
tions due to Dokshitzer , Gribov and Lipatov ^(| and 
Altarelli and Parisi [l7| are modified. The MRST group 
obtains 
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In Eq. the right hand side of the schematic evolution 
equations represents a convolution of the splitting func- 
tions with the quark and gluon distributions (which have 
an explicit dependence on the factorization scale param- 
eter /i 2 ). Inclusion of the electromagnetic contribution 
to QCD evolution introduces a "photon parton distribu- 
tion" 7(x,/i 2 ) which is coupled to the quark and gluon 
distributions. The new splitting functions that occur in 
Eq.Elare related to the standard QCD splitting functions 
by 
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Conservation of momentum is assured by the relation 



dx x 



^2 Qi(x, y 2 ) + g(x, /i 2 ) + 7(sc, y 2 



1 (5) 



It is necessary to simplify Eqs.|3 First, since the EM 
interaction is not asymptotically free, it is not clear how 
to set the starting values for the various PDFs that are 
coupled by these QED effects. In particular, it is not 
clear where the QED effects should be assumed to vanish. 
Second, inclusion of the QED couplings could in princi- 
ple more than double the number of parton distribution 
functions (one must now differentiate between proton and 
neutron PDFs, in addition to the new photon parton 
distributions). Two groups have adopted somewhat dif- 
ferent strategies, with similar overall results. Glueck et 



al. |I4| point out that the photon parton distribution is 
already of order a, as is clear by inspection of Eq.[3] Con- 
sequently to leading order in a they drop terms involving 
7(2;, /j, 2 ) from the right-hand side of Eq. [31 When they 
adopt the standard convention for DIS reactions of set- 
ting the scale /x 2 = Q 2 , Glueck et al. then obtain convolu- 
tion equations for the charge symmetry violating valence 
quark distributions arising from QED coupling, 
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Similar relations hold for the antiquark distributions. 
Glueck et al. assume that the average current quark mass 
TfTq , taken as 10 MeV, is the kinematical lower bound for 
a photon emitted by a quark. This is analogous to tak- 
ing the electron mass as the lower limit for radiation of 
photons in the earliest calculations of the Lamb shift (be- 
fore the advent of renormalization group arguments) |18| . 
Eq. [HI is then integrated from m 2 to Q 2 . The rationale 
here is to evaluate QED evolution effects while keeping 
the QCD effects fixed. Thus, the quark distributions ap- 
pearing on the right hand side of Eq. [S[ are taken from 
the GRV leading-order parton distributions [l9|. In the 
resulting integrals, in the region q 2 < [i 2 LO = 0.26 GeV 2 
corresponding to momentum transfers below the input 
scale for GRV, the PDFs are "frozen," i.e., in this region 
they are assumed to be equal to their value at the input 
scale [i% - 

The resulting valence isospin asymmetries xSu v and 
x5d v are plotted in Fig.^at Q 2 = 10 GeV 2 . For compari- 
son, they are plotted along with the valence quark isospin 
asymmetries obtained by Rodionov et al. 0, |2(j . The 
latter CSV distributions were obtained from bag model 
calculations, where charge symmetry violation was as- 
sumed to arise from mass differences of the residual di- 
quarks Sfh = nidd — rn uu and from the target nucleon 
mass difference SM = M n — M p . The quantity Sfh was 
taken as 4 MeV 0, ^J, which includes an estimate of 
the EM contribution to this mass difference. While the 
quantity Su v is quite similar in both sign and magnitude 
for both the bag model and the QED calculations, the 
QED results for Sd v are roughly half as large as the bag 
model results. As noted previously, the bag model re- 
sults for valence quark CSV are extremely close to those 
obtained by MRST using the phenomenological form of 
Eq. for the best-fit value k = —0.2. 

The MRST group solves the evolution equations 
of Eq. [31 with assumptions about the parton distributions 
at the starting scale Qq = 1 GeV 2 . At the starting scale, 
the sea quark and gluon distributions are assumed to be 
isospin symmetric. The starting photon parton distri- 
butions are taken as those due to one-photon radiation 
from valence quarks in leading-logarithm approximation, 
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FIG. 1: The isospin-violating majority xSu v (solid curve) and 
minority xSd v (dashed cu rve) valence parton distributions 
obtained by Glueck et al. 14] at Q 2 = 10 GeV 2 , assuming 
QED evolution from a scale set by the current quark mass. 
These are compared with majority (solid points) and minor- 
ity (open circles) CSV distributions obtained from theoretical 
bag model calculations [To| . 



evolved from current quark masses m u = 6 MeV and 
raid = 10 MeV to Qq. This produces different photon 
PDFs for neutron and proton at the starting scale. En- 
forcing overall quark momentum conservation from Eq.[S] 
requires valence quark isospin asymmetry at the starting 
scale. MRST assume that this takes the form 



2« - <) = e« ~ 2<) 



(7) 



Eq. is a simple phcnomcnological form chosen to obey 
the valence quark normalization condition, and it pro- 
duces isospin violating distributions that resemble the 
valence PDFs at large and small x. The parameter e is 
determined from the overall quark momentum conserva- 
tion condition. 

Having determined the starting distributions for the 
photon parton distribution, and the asymmetry param- 
eter e, the proton's quark and gluon distributions at the 
starting scale Qq are determined in the same way as for 
other MRST global fits. The only change is that separate 
DGLAP evolution equations are used for partons in the 
neutron and proton. 

If we adopt the MRST functional form for charge sym- 
metry violating PDFs, we can estimate the magnitude of 
effects one might expect in a dedicated experiment de- 
signed to test parton charge symmetry. In a recent paper 
|2l| , we estimated the magnitude of effects in two promis- 
ing experiments. The first was a comparison of Drell-Yan 
cross sections induced by charged pions on an isoscalar 



target (e.g., the deuteron). The second experiment in- 
volved semi-inclusive deep inelastic scattering involving 
charged pion production in e — D interactions. 

In this report, we consider another possible experimen- 
tal test of parton charge symmetry. This involves mea- 
surements of W production at hadron colliders, specifi- 
cally VF-boson production in high energy p— D collisions. 
This was initially suggested by Vigdor 22]. Boros et 
al. made estimates of the effects that might be expected 
at colliders such as RHIC and LHC |2^], and concluded 
that one might expect several percent effects in certain 
observables. However, these effects occurred because the 
authors had assumed very large charge symmetry vio- 
lation in the parton sea. This large sea quark CSV was 
necessary to account for significant discrepancies between 
the F 2 structure functions extracted from high energy 
U — D interactions measured by the NMC Collaboration 
and the F 2 from v-Fe DIS measured by CCFR 
26] . However, these discrepancies di sappea red when the 
neutrino reactions were re-analyzed [271 128| . 

There were two reasons for significant changes of F% 
upon re-analysis. First, experimental neutrino cross- 
sections measure a combination of F 2 and xFj,. In the 
initial analysis the quantity xF^ was calculated from phe- 
nomenological PDFs. In the re-analysis, this quantity 
was extracted from experiment, by using the fact that 
the two structure functions have different y dependences 
in the cross sections. The value of xF^ that was ex- 
tracted differed considerably from the phenomenologi- 
cal xF$, and this subsequently changed the value of F2 
that was extracted. The second significant change arose 
through the use of next-to-lea ding order (NLO) equations 
for charm quark mass effects |28|. rather than the "slow 
re-scaling" model [HHO]- Although the MRST analysis 
found evidence for sea quark CSV 0, it was consider- 
ably smaller than that extracted from the original data, 
including the slow re-scaling contribution. 

The cross sections for the processes p + D — > W + + X 
and p + D — > W~ + X have the form 



da 

dx F 



{ P D -» W+X) ~ cos 2 ® c {u{x 1 ){u{x 2 ) + d{x 2 ) 



-Su(x 2 )) + d(xi) (u(x 2 ) + d(x 2 ) - 8d(x 2 ))] + 

sin 2 C [2u(xi)s(x 2 ) + s(xi) (u(x 2 ) + d(x 2 ) - 5d(x 2 ))} 



da 

dx F 



(pD -» W~X) - cos 2 O c [u( Xl )(u(x 2 ) + d{x 2 ) 



—5u(x 2 )) + d(xi) (u(x 2 ) + d(x 2 ) — 8d(x 2 ))] + 
sin 2 6 C \2u(xi)s(x 2 ) + s(xi) (u(x 2 ) + d(x 2 ) - Sd(x 2 ))] 

(8) 

In the absence of CSV terms, if we take the sum of the 
W + and W~ cross sections, 



as(x F ) 



da 

dx P 
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(9) 



then the Cabibbo favored terms in a s are invariant under 
the exchange x\ <-> x 2l or alternatively under the trans- 
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formation x F — > —x F , where x F = xi—X2- Consequently, 
we define the forward-backward asymmetry A(x F ) as 



A(x F ) = 



a s (x F ) - a s (-x F ) 
a s {x F ) + a s (-x F ) 



(10) 



The only terms remaining in the quantity A(x F ) are 
charge symmetry violating terms, plus terms containing 
strange quarks in the Cabibbo-unfavored sector. 

We have calculated the effects to be expected for the 
forward-backward asymmetry A{x F ) for VK-production 
at a hadron collider, using the charge symmetry violat- 
ing PDFs calculated by the MRST group. We have used 
PDFs corresponding to three different sources of charge 
symmetry violation. First, we used the valence quark 
and sea quark CSV PDFs extracted bythe MRST group 
from global fits to high energy data Q. Then we have 
added the CSV PDFs calculated by the MRST group by 
including the "QED" contributions to QCD evolution, 
with assumptions about the parton distributions at the 
starting scale Q 2 = 1 GeV 2 [H ^ ■ Note that the vari- 
ous CSV PDFs were extracted using different procedures. 
The valence and sea quark PDFs were calculated in inde- 
pendent global fits to high energy data. In these fits, one 
assumed a particular functional form for the valence (or 
sea) quark CSV PDFs, which depended upon an over- 
all variable parameter for the strength. That parameter 
was determined by minimizing the \ 2 °f the global fit. 
For simplicity MRST neglected the Q 2 dependence of the 
CSV distributions. 

In these various global fits, one obtains different va- 
lence parton distributions in the minimization process 
(the sea quark and gluon distributions were essentially 
identical to those obtained assuming charge symmetry). 
So the best fit valence quark PDFs obtained by MRST 
when they allowed valence quark CSV differed somewhat 
from those obtained when they allowed sea quark CSV. 
In addition, the MRST global fits that allowed parton 
CSV did not explicitly include the "QED" contributions 
to QCD evolution. We explore the various contributions 
to the W production asymmetry, despite some questions 
regarding the consistency in the different parton distri- 
butions that give rise to CSV effects. 

In Fig. [21 we plot the forward-backward asymme- 
try expected for W production, as defined in Eq. EH 
The parton distribution functions are obtained from the 
MRST analysis that includes the "QED" contribution 
to DGLAP evolution; this analysis includes electromag- 
netic couplings in the evolution equations that give rise 
to isospin violation p| , as given by Eq. |3J The input 
data for these global fits was that used in the MRST2004 
analysis [13 ■ The top figure is calculated for y/s = 500 
GeV, and the bottom figure is calculated for ^/s = 1000 
GeV. All three curves in Fig.[21mclude the sea quark CSV 
terms from Eq. [21 They differ in the amount of valence 
quark CSV of the type given by Eq. ^ I n eacn graph, 
the solid line corresponds to no valence CSV (k = 0), 
while the short-dashed and long-dashed curves are for 
k = —0.8 and n — +0.65, the parameters corresponding 
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FIG. 2: The forward-backward asymmetry A(x F ) defined in 
Eq. il l Oil as a function of x F . Top graph: ^fs = 500 GeV; 
bottom graph: y/s = 1000 GeV. The curves include CSV 
generated by QED effects, and sea quark CSV described by 
Eq.|U They differ in the amount of valence quark CSV defined 
by Eq. Solid curve: no valence quark CSV, ft = 0; long- 
dashed curve: k = +0.65; short-dashed curve: k — —0.8. 



to the 90% confidence limit obtained in the MRST global 
fit. 

Fig. |21 shows that the sea quark CSV and the QED 
isospin violation tend to produce negative forward- 
backward W asymmetries. The contributions from va- 
lence quark CSV, with magnitudes at the 90% confi- 
dence level extracted by MRST, contribute a roughly 
equal magnitude to the asymmetries produced by the 
other sources of CSV. For negative values of k, which 
agree with theoretical estimates of valence quark CSV 
the valence CSV terms tend to cancel the asym- 
metry produced by sea quark CSV; while for positive 
values of n the various sources of isospin violation tend 
to add. Note that the predicted forward-backward asym- 
metries are quite small. The magnitude is less than 0.01 
for almost all values of x F . These results are consider- 
ably smaller than those obtained by Boros et al. [23|. 
who predicted rather large positive values for A(x F ), as 
large as A(x F ) ~ +0.07 for x F ~ 0.7. There are two rea- 
sons for this difference. First, the sea quark CSV terms 
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are substantially smaller for MRST than those extracted 
by Boros et al, by a factor of five or six; the sea quark 
CSV terms obtained by MRST and Boros also have op- 
posite signs. Second, with the very large values of sea 
quark CSV extracted by Boros, the Cabibbo unfavored 
terms were negligible. However, with the much smaller 
sea quark CSV obtained by MRST, the Cabibbo unfa- 
vored terms can no longer be neglected, and they tend 
to cancel the Cabibbo favored contribution. 

This calculation of the forward-backward asymmetry 
for W production finds that the expected effects are quite 
small, generally less than one percent for the range of 
CSV effects obtained by MRST. With the expected lev- 



els of isospin violation in PDFs, it will be necessary to 
measure these asymmetries to better than one percent if 
this observable is to provide a test for charge symmetry 
violation in parton distributions. 
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